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The problem of the generation of an electromagnetic pulse by an expanding, infinitely 
conducting, spherical plasma under the earth is considered. The solution consists of the 
derivation of an appropriate tensor Green's function for the half-space which reduces the 
problem of determining the electric field at any point in space due to the current density 
generated by the interaction of the plasma with the earth's static magnetic field to evaluating 
an integral. The vertical component of the electric field at the earth's surface which is 
generated by the mechanism is calculated. 



1. Introduction 

The problems of determining the electromagnetic field resulting from a vertical and a 
horizontal dipole imbedded in a conducting half-space have been solved in the past essentially 
by following Sommerfeld's method for solving half-space boundary value problems. A similar 
problem such as the one of determining the field due to an arbitrary vector current distribution 
J(r, t) imbedded in a conducting half-space with a time independent, uniform conductivity a, 
is often of some interest. Such a problem, for example, may arise from the electromagnetic 
radiation of an expanding plasma which may result from a very intense underground explosion 
[Taylor, 1950]. The electromagnetic field generated by the current density J can be obtained 
by first evaluating the appropriate Green's function which satisfies the boundary conditions 
and then evaluating the volume integral of the product of J and the Green's function. In 
this paper, the Green's function for a time-varying source is derived by making use of Fourier 
transforms. The field expression at any point in space is then obtained. Numerical results are 
given for the vertical component of the electric field at the earth's surface. 

2. Formulation of the Problem 

From Maxwell's equations it is seen that the electric field E(r y co) as a function of position r 
and angular frequency co satisfies the equation 

VXVX#(?, co)-/3 2 S(7, co)=ico Mo 7(7, co) (1) 

where E(j, co) and J(,7 co) are respectively the Fourier transforms of the corresponding time- 
varying electric field E(r, t) and the current J(r, t) , denned by 

E(r,«)=4= P E(r, t)e«"dt, J(r, «)=-= P J(r, t)e*»'dt 

■yjZW J -oo -y/Zir J -oo 



i The work reported here was partially sponsored by the Air Force Systems Command, Rome Air Development Center .Griffiss Air Force Ba 
\rw York, and by the Advanced Research Projects Agency, Washington, D.C., 20301, under Contract AF 30(602)-2500. 
2 Now ;it American Nucleonics Corporation, Glendale, Calif. 

147 



and p is the propagation constant defined by 

2 f/x € co 2 in the nonconducting half-space 
\ju e co 2 -\-iuixa in the conducting half-space. 

The symbols /x , to, and a denote respectively the free-space permeability, free-space dielectric 
constant and the conductivity of the conducting half-space. 

Thus, if a dyad r is defined which satisfies the equation 

VXVXf-i3 2 f=<5(7-7 / )/ (3) 

where 7 is a unit dyad and V satisfies the appropriate boundary conditions at the interface of the 
finitely conducting and nonconducting half-spaces, the electric field intensity can be written as 

E=i(afjL Q | Y-Jdv', (4) 

the integration being taken over the region of space V where J does not vanish identically. 

Let the conducting and nonconducting parts of the space media be respectively referred 
to as 1 and 2, as shown schematically in figure 1. Also let a cartesian coordinate system be 
chosen such that the z axis is normal to the interface between the media 1 and 2. 

Since (3) has coefficients which do not depend on x and y, one may introduce the Fourier 
transform 



rfe Vj * )= Gb) / J— T(x ' y> z)enx+i " vdxdy 



(5) 



into this equation. Then the operator V is transformed into V^ where 



V F =—i&—ivy+fa z- (6) 

In (6) and in the following, a caret is used to indicate a unit vector in the respective coordinate 
direction. Thus (3) can be written as 

V F XV F xF(£, rj, 0)-/3 2 r(£, r,, z)=l±8(z-z')e i t x '+W. (7) 

Equation (7) represents a system of ordinary differential equations, in the variable z, for the 
components of I\ Appropriate boundary conditions must be satisfied. The continuity of the 
tangential components of K for an arbitrary current source means that r rr , T TZ , and Tee are 
continuous at 2=0. Continuity of the tangential component of H (proportional to V X f) 
implies that the rO, Or , and Bz components of V X r are continuous at 2=0. The other tan- 
gential components of V X r vanish identically. 



MEDIUM (I) 



• r' = (x', y' , z' ) 

Figure 1. Coordinate system used in deriving 
Green's function. 
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The solution of (7) satisfying these conditions is readily obtained. If the source point is 

in medium 1 and the field point in medium 2, the Green's function r which satisfies (3) and 
this boundary condition can now be expressed as 
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(8) 



where 



and the subscripts on the propagation constants refer to the conducting and nonconducting 
media. 

3. Pulse Generated by an Expanding Plasma 

This Green's function will now be used to solve the problem of the emission of an electro- 
magnetic pulse by a very intense explosion under the earth. It will be assumed that the 
explosion creates an infinitely conducting, expanding spherical plasma which interacts with the 
earth's magnetic field. It is known that the magnetic field is frozen in the plasma as it ex- 
pands [Cowling, 1957]. Thus, if the expansion is assumed to begin from a point, no magnetic 
flux will exist inside the volume containing the plasma. 

The physical mechanism for canceling the earth's magnetic field, J1 Q , inside the plasma is 
a surface current which is induced by B . Hence, the plasma may be replaced, for computa- 
tional purposes, by an equivalent current density J in the earth. This is more convenient 
than solving directly the difficult moving boundary value problem which requires, in addition 
to the vanishing of the normal component of the total magnetic field at the plasma surface, 
the vanishing of the tangential component of 3 

at the plasma surface where UJ=E(r, t) and B=B(r, t) are respectively the electric and magnetic 
fields generated by the plasma. Here a=a(t)r is the radius of the plasma. 

To obtain an expression for the equivalent current density generated by the expanding 
plasma, we first remark that if the explosion takes place at a sufficiently great depth, the effect 
of the earth-air boundary on the field near the plasma will be negligible. In fact, even ne- 
glecting attenuation in propagation through the earth and assuming that the earth-air boundary 
is a perfectly reflecting plane, the ratio of the reflected field to the primary field at the surface 



3 To see that this boundary condition must be satisfied instead of the usual vanishing of the tangential component of E, one can consider, in 
addition to the fixed coordinate system used above, a coordinate system which momentarily is moving with the velocity of some point on the 
surface of the sphere. In this coordinate system, a portion of the plasma surface is stationary and the boundary condition is that the tangential 
component of E vanishes. Transforming back to the fixed coordinate system, we find that the boundary condition is that stated above. 
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of the plasma is of the order of 

2d ) 

for low frequency waves where d is the depth of explosion origin and B max is the maximum 
radius which the plasma achieves. For example, the effect of the earth-air boundary is less 
than 0.1 percent for a plasma which achieves a final radius of 25 meters if d=150 meters. 

Thus, to calculate the current density at the surface of the plasma which will cancel the 
earth's magnetic field inside, it can be assumed that the plasma is located in an infinitely 
extended medium with conductivity a. Let r, d f , (/>' be the spherical coordinates of the point 
r in a coordinate system using B as the polar axis and the center of the plasma as the origin. 

In this case, the symmetry of the problem dictates that the current density be in the 4>' direc- 
tion and proportional to sin 6'. For, it is clear that any vector associated with the electro- 
magnetic field is independent of the angle </>'. Furthermore, if a reflection in any 
plane containing B followed by multiplication by minus one is performed, it is seen that 
B Q (which transforms as a pseudo vector) and the geometry remain invariant. Therefore, all 
vectors associated with the EM field must be invariant under these operations. This implies 

A A A 

that vectors can have only a <j>' component and pseudovectors r and 6' components. The 
fact that the angular dependence of J is given by sin 6' can be seen by recalling that a plasma 
moving in a radial direction interacts with the external field B by a term proportional to 
v X Bo (where v is the radial velocity vector) . By substituting a driving term of this form into 
Maxwell's equations and using the above results, it is immediately found that J is propor- 
tional to sin Q f '. 

Hence, J can be written as 

J=K(t) sin d'8(r—a)4' (9) 

where the 5 function has been inserted because the current is confined to the surface of the 
plasma. Here a=a(t) is the radius of the plasma and K(t), which may depend on time, is 
to be found. 

To determine K, we find the magnetic field which is generated in a medium of conductivity 
<r by a current density of the form indicated in (9) . It is convenient to Fourier transform into 
the frequency domain. Since the electric field Green's function in an infinite medium can be 
expressed as [Morse and Feshbach, 1953a] 

t= W~l ^0l~-~'l oo i o/_i_i (7_™v . 

I+ W] NT* -S.S S WW) e » h+tff • {^^ > )M< ?H (^ < )+iV-(^>)iV-(/3r < )} 

even 

where M and N are vector spherical harmonics, e :n =< ., ~ and the symbols r> and r < 

indicate respectively the larger and smaller of r, r ', it is found, as a result of the orthogonality 
relations among the vector spherical harmonics, that 

B(r, o>)=±VXE=i^ 2 Nlll Q<h0 (M rih(Pr')J(r', u)r'Hr' +N<*L. li0 (' ji(Pr')J(r' 9 u)r n dr f 

^ Jr JO 

where h x andji are spherical Bessel functions and J(v' ', co) is the Fourier transform of K(t)8(r—a). 
Fourier transforming to the time domain, it is found, after taking (9) into account, that 



z«e-' w( '-''>i/xo|8 2 a 2 



{Mil, i, oiMhi (pa)H(a-r) + iV«> en , ,, „(/3r) j, (fia)H(r-a)} (10) 
where the argument of a is t' . 
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Now, it is readily seen from dimensional arguments [Cowling, 1957] that if a non-uniform 
spatial distribution of the magnetic field exists inside the sphere, the time for the non-uniformity 
to decay to zero is of the order of /x ot 2 . For a sphere with a radius of 25 meters and an earth 
conductivity of 10~ 2 mho/m, this time is about 10" 5 seconds which is much smaller than any 
practical expansion time for the plasma. Thus, the field at any point within the radius a(i) 
can be considered to be quasi-static and equal to the field at r=0. Hence, to evaluate K(f), 
we may set r=0 in the following. 

In this case, (10) reduces to 



5(0, t)=± J" ' dt'K(t')iijtfa*N£L. i.o(0)Ai(j8a) 



~2tt 3 



\ (cos B'r— sin 6'$') f d w «- , » (| - | ' ) [/So+i]e- Ma - (11) 

*~* %J — 00 



Since /3a is small for all frequencies of interest, the term e i?a may be expanded in a power series. 
Keeping terms of order j3 2 a 2 , it is found that 

/7«), *)=g| (coBtf'r-Bin 9'$') C dt'K{t') f " **«-'"«-"> fl+^1 

= Mo | (cos 0'r-sin 6'$') f'dt'K(t') [*(«-*')+^ { ~W ^-*>* Jp} «(*-*')] 

=Mo | (cos fl'r-sh, «'*') [tf(0+*f ^ (#« 2 )-f £ (Xo 1 )} 

It is seen, then, that the magnetic field generated by the current density e/ will cancel the 
earth's magnetic field 

£ (cos B'r— sin 6'6 f ) 
if K(t) satisfies the differential equation 

-§ ^=K(t)Vf j t (Ka>)-f Jp (^ 2 )- (12) 

In particular, for slowly expanding plasmas, K(t) is nearly a constant: 

^ Mo 

j 2 

This will be true in cases of practical interest, for the term -~ -j- (Ka 2 ) is of order ° smaller 

than the leading term in (12). This is very small. The third term of (12) is of even smaller 
order. 

Thus, the current density 

J(r, 0=-i — sin B'b(r'-a)4>' (13) 

2 Mo 

will approximately cancel the earth's magnetic field inside the plasma. It can be readily 
checked by expanding the integrand in (10) in powers of fir that this expression for the current 
density will lead to a field B(a, t) at the boundary of the plasma whose normal component differs 
from — cosd'Bo by a term in the integrand of order (13a) 2 smaller than the leading term. A 
similar calculation shows that the second boundary condition is also satisfied to the same order 
of magnitude. For frequencies near 1 kc/s, an earth conductivity of 10" 2 mho/m, and a maxi- 
mum sphere radius of 25 meters, the error of order (fia) 2 in the integrand is about three percent 
of the leading term. 
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MEDIUM (2) 




EXPANDING 
PLASMA 



Figure 2. Illustrating geometry for electric field 
calculation for an expanding plasma. 



To calculate the electric field at the earth's surface produced by the expanding plasma, 
it is necessary to write (13) in terms of the coordinate system in which the Green's tensor is 
expressed. Let the earth's magnetic field be inclined at an angle 12 to the vertical axis, as is 
illustrated in figure 2. The primed system of coordinates is also illustrated with the x and x' 
axes oriented so that they coincide. Then the y' axis is in the y— z plane. The following- 
equations hold: 

A/ A 

X X 

y' = cos Six— sin 12 z. 

A 

Furthermore, the vector </>' which occurs in the expression for the current density is given by 

#' = — sin ^/z' + cos <j>'y' — — sin 0'x + cos <// cos 12 y — cos 0' sin 12z. 
We want to express the quantity sin 6'$' in terms of the unprimed variables. Now, 
sin 0' sin <!>'=£' *y=sin 6 sin </> cos 12— cos 6 sin 12 



and 
Thus, 



sin 0' cos <// = *'.,£' = sin cos <j>. 



sin 0'$>' = [— sin 6 sin </> cos 12+cos sin £2]x+sin $ cos cj> cos 12?/ — sin cos </> sin 122. (14) 

The vertical component of the electric field at the earth's surface can be found from (4), 
(8), (13), and (14). It is convenient to perform jthe volume integral indicated in (4) before 
attempting the £, rj integrals in the expression for T. Thus, we write 



where 



Ez= ~±?f\ r. ^ft+^^-^+^bj-^^^-^^ii (is) 

p«= f J«(co)e^+^ +7 i< 2+d W. (16) 

J sphere 



In this equation a=x, y, or z and J«(co) is the Fourier transform of J a (t) . We have also set z=0 
in the equation for E z and have assumed that the center of the expanding plasma is at a depth 
d below the interface. 
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The integrals in (16) are to be evaluated using (13) and (14). Then p Z} for example, is 
given by 

p x =~ — f F[5(r-a)]r 2 dr f * dd sin 6e y >- r "» e 

* Mo J Jo 

'27T 

d<t>[— sin sin cos 12+cos sin a], e <x' ■**««>« (*-<*> 



•/; 



where jP[5(r— a)] is the time Fourier transform of 8(r—a(t)) and we have introduced the spher- 
ical coordinates r, 0, <j> centered at the plasma center {z axis vertical) . The parameters X and a 
are defined by the equations 

£=X cos a 

r]=\ sin a. 

The <j) integration is readily performed and we obtain 

Vx=~\— ( ' F[8(r-a)]r*dr (' dde^™ 9 
2 Mo J Jo 

•[— sin 2 cos 1227ri sin aJi(\r sin 0)-fsin cos sin tt2wj (\ sin 0)] 

where J and Ji denote Bessel functions of the zeroth and first order. To perform the inte- 
gration, we use the following result listed in Morse and Feshbach [19531)]: 

* e iz cob6 cob uj m ( z sill sin ^)pj*( cos u) sin udu=i n - m J—P%(cos B)J n+m {z) 

where P™ is an associated Legendre polynomial. By using this result, one finds for p x% and in 
the same manner, p y and p 2 

p x =-~ [-iXsinacosQ+TxsinO] |^2 f^^r-a)]^/^)^ 

Pi - Mo J 

p y = — ~ cos a cos 12 - — I F[8(r— a)]ji(fiir)r*dr 

Pi * Mo J 

p z =^ sin 1] cos a I -° f ^(r-aJtofor^dr 

Pi * MO J 

where ^ is the spherical Bessel function of order one. 

Substituting these results into (15), we find, after performing the a integration and making 

the substitutions p 2 =x 2 -\-y 2 , tan B=- and using the fact that ^— Jo(\p) = —Ji(\p)\ that 

S z (p,^rf, w ) = -2i«/io|8isinacos^|^. f^Cr-aJb'iOSxrJr 2 ^!- f" * *'* «7 (Xp)XriX (17) 

^ Mo J Op Jo 7iP2+72Pi 

The integral over X in (17) represents the effect of propagation over the earth and has 
occurred in many studies of radio propagation over a conducting half-space. For distances 

p<50 km (but much larger than the skin depth Im /Sr 1 ) and for frequencies /=^-<10 3 cps, 

AW 

one can show that [Space-General Corporation, 1960] 

?) r°° p-yi* piPid i- -in pifap 

V" -WT—iv e/o(Xp)X^X « V -fa+~ — (18) 

dp Jo 7i$+Y 2 /3i ffl L PJ P 

correct to within 15 percent. 
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To evaluate the remaining integral in (17), it is necessary to be more specific about the 
nature of the source. According to Taylor [1950], the radius of the plasma increases with 
time (at least during the initial instants) according to the formula 



a(t)=At 2/s (0<t) 



(19) 
x at 



where A is a constant. We will assume that this expansion continues until a time t } 
which time the plasma radius is R and then stops abruptly. 

With this expression for a(t), it is convenient to perform the r integration before per- 
forming the Fourier transform in the source integral in (17). Thus, 



Jo y27r Jo Jo 



drr'jxiMKr-ait)) 



tfWjxfaaWe^dt 



_ i r 

■yj%r Jo 



The last form results as a consequence of our assumption that a(t) =R for t>t mSLX . The quan- 
tity e is a small positive number whose purpose is to shift the pole off of the axis of integration 

when . is multiplied by a function and integrated over co. Because of the multiplying 

factor co in (17), no improper integrals will occur in our problem and we can set e=0 in all that 
follows. 

Because of the attenuation factor e~ Im ^ d which occurs in (18), high frequencies will be 
attenuated very rapidly in the earth. Furthermore, the maximum radius R of the plasma 




IO" 3 
TIME (SECONDS) 

Figure 3. Electric field strength vs time after be- 
ginning of expansion. 
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where 



will not be large in most conceivable cases. Thus fti^l and we can expand the spherical 
Bessel functions ji(Pid) and ^ (ft 72) in a power series, keeping only the first term. When this 
is done, one can write, after collecting the above results 

i p -fop r i ~i 

E z (p, 6, d, u)=-== B R* sin Qe*W ~ift+- A(<a) cos 6 

-y 2tt P I— P-l 

Yl(co)=ie i " ! ^+co ( lm * X (j>) e tat dt. (20) 

Finally, the time dependent field is 

E Z { P , e, d, 0=4= r E,( P , e, d, uie-^da. (21) 

■\2iTT J -co 

Figure 3 illustrates some results obtained by use of (20) and (21). The conductivity of 
the earth was set to cr = 10~ 2 mho/m. B was taken equal to 0.5 X10~ 4 webers/m 2 and sin 
12 = 0.707. The depth of the plasma was taken as <i= 165 meters and the plasma was assumed 
to expand to a radius R=26 meters in a time £ max = 1.8X 10~ 3 sec while obeying the law expressed 
by (19). The time integral was evaluated numerically on an IBM 704 computer. Figure 3 
presents the vertical component of the electric field pulse at the earth's surface for several 
distances p evaluated at 0=0. 

4. Conclusion 

The tensor Green's function for the electric field in a conducting half-space has been ob- 
bained. With the aid of this function the calculation of the electric field generated by an 
arbitrary current distribution imbedded in a conducting half-space is reduced to evaluating an 
integral. The formalism was illustrated by an example which indicated the usefulness of this 
approach. 
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